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Recently several numerical results have been reported in the three and four- 
dimensional simplicial quantum gravities concerning their phase structure [Q, 0. Espe- 
cially in four dimensions, it is suggested in |0] that there is a second order phase transition 
point indicating the continuum limit of the model. In view of the importance of this 
suggestion, it is desirable to develop an alternative analytical approach to the dynamical 
triangulations. 

In two dimensions, the matrix model has been found a useful tool in studying the 
simplicial gravity, and there have been some attempts to extend it to higher dimensions . 
In this so-called tensor model, one considers an integral over rank-3 tensors rather than 
matrices, and the perturbative expansion of the integral generates randomly triangulated 
three-dimensional manifolds. However, in this model, it is found difficult to have control 
over topologies generated by the random triangulation. This is because the model does 
not contain the sufficient number of parameters to distinguish different topologies. To 
solve this problem, Boulatov []5| proposed an improved version of the model defined in 
terms of functions with three variables on a group or a quantum group rather than the 
tensors. In his model, the perturbative expansion of the integral is expressed as a sum over 
triangulated three-manifolds weighted with some topological invariant which depends on 
the choice of the group G. There is a large variety of topological weights corresponding to 
various choices of groups, and one may hope to gain control in summing over topologies 
with suitable choice of groups. 

The purpose of this paper is two-folded. First we will extend the Boulatov's model 
to four dimensions and check that its perturbative expansion is expressed as a sum over 
triangulated four-manifolds. We then show that the weight for the summation is invariant 
under the Alexander moves of the simplicial complex, and thus depends only on 
the combinatorially equivalent class of the triangulation. This model could be useful for 
analytical study of the four-dimensional simplicial gravity. It is known that there is no 
algorithmical way to classify topologies in four dimensions, but it would still be possible 
that, for a suitable choice of G, the perturbative expansion of the model is dominated by 
S'^ or by some restricted set of topologies i. 

^ This seems to be the case in three dimensions [^. 
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The topological invariants obtained in this way are interesting in their own right, and 
examining their properties is the second purpose of this paper S. In order to solve the 
quantum gravity in four dimensions, independently of how one defines it, one would need 
to have a good understanding on four-dimensional geometries. Thus it should be useful to 
study various geometrical invariants formulated in the language of quantum field theory 
and statistical mechanics. In the three-dimensional model of Boulatov, if we choose G to be 
the quantum group SUq{2), the weight in summing over topologies is constructed from the 
6j-symbols of SUq{2) P]. The weight then reproduces the three-dimensional topological 
invariant defined by Turaev and Viro |T^, whose q ^ 1 limit has been studied many 
years ago by Ponzano and Regge It has been shown in |]12| that the Turaev-Viro 
invariant is related to the partition function of the Chern-Simons gauge theory [jl5 
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three dimensions. We shall see that this feature persists in four dimensions. For q = 1, 
the four-dimensional model constructed here makes use of the 15j-symbols as well as the 
6j-symbols, and it may be regarded as a four- dimensional version of the Ponzano-Regge 
model. It is also possible define the g-analogue of this model, whose partition function 
could be regarded as a four-dimensional version of the Turaev-Viro invariant. Although 
the proof of the topological invariance given here only applies to the case of = 1, we 
present an evidence that it should also hold for when g is a root of unity. For q = 1, the 
lattice model is related to the so-called BF model [|16[ [|l3 We have not been able 

to identify a quantum field theory corresponding to the g-analogue model, and this task is 
left for future investigation. 

In the three-dimensional Chern-Simons theory, the Wilson operators along closed 
loops make a basic set of observables and their expectation values are related to the Jones 



invariants of knots [|I4[ . One may also consider the 't Hooft operator , but it is known 
that the 't Hooft operator is identical to the Wilson operator in the Chern-Simons theory 
|21[j . In four dimensions, on the other hand, it is a surface that makes a non-trivial knot. 
We study the 't Hooft operator in the BF model, which is associated to S"^ embedded in 
four dimensions and differs from the Wilson operator, and discuss how to realize it in the 
lattice model. Both the 't Hooft operator and the Wilson operator create physical states 



^ In this paper, we discuss geometries in the piecewise linear category. If two differential 
manifolds are diffeomorphic, their smooth triangulations are combinatorially equivalent Q . Thus 
an invariant of the piecewise linear topology gives a diffeomorphism invariant of smooth manifolds. 
In four dimensions, it is not know whether a homeomorphism of two combinatorial manifolds also 
implies their combinatorial equivalence. This is known to be the case in two and three dimensions. 
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on S'^ X embedded in the four-dimensional manifold M. In the case of the 't Hooft 
operator, S'^ x 5"^ is a boundary of 5"^ x embedded in M, and its S'^-part is contractable. 
On the other hand, for the Wilson operator, S'^ x S"^ is a boundary of x S'^, and the 
5'^-part may link with other operators. Correlation functions of the these operators in the 
g-deformed version of our lattice model would define a new type of invariants of knots and 
links in four dimensions. 

Let us define our model. We start with a real-valued function of four variables 
4>{.dii 92-, g^, di) on G {gi G G). For simplicity of notations, we describe the case when 
G = SU{2), but extensions to other groups should be obvious. If we take G to be a finite 
group, the model constructed below may be regarded as a four-dimensional version of the 
Dijkgraaf-Witten model . The function (p can be expanded in terms of matrix elements 
^mnig) °f ^ spin-j representation as 

Hgi^ g2. gs. g4) = Yl '^m;n'...m4n4^mini(fi) ■ ■ ■^^4^4(^4)- (1) 

ji,m,i,ni 

We then require (p to be invariant under the right action of G, 

(pigiU, g2U, gsU, g^U) = (j){gi, g2, 5^3, fi'4) , {U e G) 

This condition is equivalent to 

4>{gi. g2, gs: gd = j dU(/}{giU,g2U,g3U,g^U), (2) 

where dU is the invariant Haar measure on G normalized as J dU = 1. If G is a finite group, 
the integral is replaced by a summation over group elements divided by the cardinality of 
the group. With this normalization, integrals of products of -D^n's become 



3 

' Ji 32 J3 \ ( Jl 32 J3 



3 



mi 7712 I \ni n2 ns 



(3) 



2=1 j ,Tn ,m,' ,n ,Ti' 



Jl J2 J ] J J3 J4 

^ \ m' 7713 1714 



Jl J2 J \ nn' (J J3 J4 

ni n2 n J ^ \ n' ns 714 



Here we used the Wigner 3j-symbol which is related to the Clebsch-Gordan coefficient as 

71 79 7-? \ (-l)ji-h-m3 

J I J I J6 \ _\ ■{jij2mim2\j3m3) 



mi 1712 rns J VSjs + 1 
and the metric defined by 

^rnra' ^ ^/2jTl 

By substituting (|1|) into (0) and using the above formulae, we see that (/> is expanded as 



(4) 



where 



mim2;m3m4 / ^ V^mi ni .. .7714714 yji yj4 X^n'i n'2 \^n' tl'^ ^4/ 



77, n' ,ni,n' 



(5) 

There are several conditions that the coefficient M should obey. The real-valuedness 
of (p implies 

E4 

""■^ 777 1777 2; 777 3 777 4 V J —7771 —7772 ; —7713 —7714 5 



where we used the hermiticity of -0^77 

^^77(^7) = (-l)^'-^^'^"'"'^^i777-77 (^7)- 

We also require 4>{gi, 92^ 93, 94) to be invariant under cyclic permutations of any three of 
its arguments. By the definition (H), this means that the coefficient satisfies 

M^^^^, = Y^i-iy^+^^+^^+^W{2j + l)(2j' + 1) j-^.J ^.J ^^^M^l^^^^^ (6) 

and similar relations obtained by cyclic permutations of (ji, 725^35 i4)- Here the symbol 
{■ ■ ■} in the right hand side is the Racah- Wigner 6j -symbol describing recombination of 
three angular momenta. 
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Now let us define our action. 

■^1=1 ' ^ 1=1 

fi'3, gs, fi'9)^(fi'9, fi'6, fi'2, fi'lo)</'(fi'lO, fi'8, fi'5, 

A four-dimensional simplex consists of five tetrahedra (3-simplexes), ten triangles (2- 
simplexes), ten links (1-simplexes) and five vertices (0-simplexes). If we imagine that 
the function (t>{gi, g2, gs, g^j is associated to a tetrahedron with the group elements on 
its faces, the kinetic term in the above action is seen as overlapping of two tetrahedra 
while the interaction term represents gluing faces of five tetrahedra to make a 4-simplex. 
Substituting (^) into the above and using the action can be expressed in terms of M 
as 



32,3334,, 3 1 2 
2 / J I — 7ftim2 ,77137714 I 

J, 771 



- A ^ (_i)E-!iO'i+"^^)(_i)ELi(^^+"0 



J, 777,77 




T\/r32l2;h',l4j3 T\j34l4,l5':l635 yi#j2 ^6 Jl y\^j3 ^8 ;Uoi5 T\/r34ll0',h;hjl 

-''^7n27l2;—7l4 7773-'^-'77l4n4;—na7775 -''^—77l2n6;—n8 777l-''^— 7713778;— nio— 777' —77147110;— 712— 7771 ' 



(7) 




^(_i)EJ!i(^'^+"^^)(_i)E!=i(^»+^ 



m,n. 



ji h k 
mi rii —77-2 



7712 n2 -ns / VTO3 ns -n^ J \m4 n4 -n^ 

is ^5 ^6 W 32 k ^7 W ji h k 
ms 77-5 -ne / \ -TO2 -ny J \-mi nr -ns 

js k \ [ -^5 k ho \ f JA ho h 

-7713 ^8 -riQ ) \-m^ ng -7710/ V -"^4 ?^io -ni 

is the 15j-symbol of the third kind ||2^]. It is not difficult to see that the 15j-symbol 
defined in this way can be expressed in terms of the 6j-symbols. Although the 15j-symbol 
does not possess the full rotational symmetry of the 4-simplex, the symmetry is recovered 
in the interaction term of the action thanks to the condition satisfied by M. 
A partition function of our model is defined as 

Z = l n dM^^l-^'.^^%,^exp{-S{M)). (8) 
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As in the case of two and three dimensions Q , we expand this integral in powers of 
A. Each term in this perturbative expansion is constructed from a propagator derived from 
the kinetic term in the action and a vertex given by the 15j-symbol. The non- vanishing 
components of the propagator are 

_ 1 4 

{Mi^^^^^ Ml',^^, ) = - n '^-^ V(2j + i)(2j' + i) 

i=l 

and others obtained by the cychc permutations of (ji, ^2, is, i4)- By combining these 
elements together, the perturbative expansion of the integral (H) is expressed as 

^ = (9) 

^ ' sym ) 

where the sum is over oriented four-dimensional combinatorial manifolds, NsymiC) is 
a rank of symmetries of C if any, A^4(C) is a number of 4-simplexes in C, and Zc is given 
by a summation which can be written symbolically as 

zc=Yi n (2j,+i) n n {i^j}. (lo) 

j t:triangles tetrahedra A— simplexes 

To describe the precise content of the summand in the above, it is useful to draw a 
closed ^"^-graph on each 4-simplex in C as follows. Each 4-simplex has a topology of a 
four-dimensional ball and its boundary is covered by five tetrahedra. On the boundary, we 
can draw a dual (^^-diagram by setting a vertex in each tetrahedron and and by connecting 
these vertices by links dual to the faces of the tetrahedra. We then split each vertex into 
two (^"^-vertices and add an extra-link between them. There are three different ways to split 
each (/)4- vertex, and we make an arbitrary choice of them at this moment. This defines 
a closed (p^-graph on each 4-simplex. By attaching a spin (j, m) on each link and the 
3j-symbol on each vertex of the (/)'^-graph S and by summing over the magnetic angular 

4 Although the 3j-symbol is the cyclic symmetric, it changes its sign under odd permutation 
of the indices. Thus, to be precise, the association of the 3j-symbol to the (j}^-vevtex is unique 
modulo the sign factor. Since each tetrahedron is shared by two 4-simplexes, this ambiguity is 
removed by choosing the same sign convention on the (^"^-graphes on the two 4-simplexes. 
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momenta m, we obtain a sum of product of the 15j-symbol and several 6j -symbols p3 



Different choice of the splittings results in different combination of the 6j-symbols. This 
procedure is applied to every 4-simplex in C in such a way that parts of (/)^-graphes of two 
neighboring 4-simplexes match on their common tetrahedron. As a result, we obtain a sum 
of product of 15j-symbols, each of which is associated to a 4-simplex in C, and 6j-symbols 
which depend on how one splits the 0^-vertices. We then perform a sum of this object 
over all the spins on the extra-links which were added when we split the ^^-vertices. By 
using the orthonormality of the 6j-symbol, one can show that the result of this summation 
is independent of the choice of splitting of the vertices and depends only on the structure 
of the complex C and the spins on the links dual to the 2-simplexes in C. This is what we 
represent by the product Hl^j} IlllSj} in the summand of (p^Ol). 

Now we would like to examine the dependence of Zq on the simplicial complex C. 
In 1930, Alexander defined a set of basic transformations of simplicial complexes and 
proved that two complexes are combinatorially equivalent if and only if they are connected 
by a sequence of these transformations. These basic transformations are called the Alexan- 
der moves. We would like to know how Zq behaves under the Alexander moves. Recently 
Gross and Varsted introduced another set of transformations which they showed to be 
equivalent to the Alexander moves for dimensions less than or equal to four 0, and we 
shall use them below. 

In four dimensions, there are five basic moves 0; 

(1) Consider a 4-simplex and add one point at its center. We then draw five 1-simplexes 
which connect this point to the five vertices of the 4-simplex. This decomposes the original 
4-simplex into five 4-simplexes. 

(2) Consider two 4-simplexes sharing one tetrahedron. There remains a pair of vertices 
which are not shared by the 4-simplexes and we connect them by a 1-simplex. This 
recombines the two 4-simplexes into three 4-simplexes. 

(3) Consider three 4-simplexes sharing one triangle, and suppose that there are three 
tetrahedron each of which is shared by two of the 4-simplexes. Then there are three 1- 
simplexes each of which belongs to only one of the 4-simplexes. The three 4-simplexes are 
recombined into three 4-simplexes sharing a triangle spanned by these three 1-simplexes. 
The other two moves are obtained by reversing the moves (1) and (2). 

To examine the behavior of Zc under these moves, it is useful to know that the 15j- 
symbols which appear in (|10|) can be expressed as a sum of product of the 6j-symbols p3[1 . 
Then the Biedenharn-EUiot identity and the orthonormality of the 6j-symbols can be used 
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to show that Zc is invariant under these moves upto a multiphcative factors. Rather than 
showing the detail of the computations, we present here another proof, similar to the one 
used in in studying the three-dimensional model. 

By using the formulae (|^) and the definition of the 15j-symbol in terms of the 3j- 
symbols, we can rewrite Zq as follows. 

Zc = Y.Y. n + n ldUTD^^^-{UT)---D^''-{UT) (11) 

j m,n t:triangles T:tetrahedra 

where each tetrahedron T carries a group element Ut, and Ji,T5 •••5 J4,T are the spins 
on the links of dual to the four faces of T. The matrix elements D^^s are multiplied 
together around triangles in C. Each triangle t is shared by a finite number of tetrahedra 
Ti^\ T2*'*, ...,Tn^ , and we can perform the sum over the spin-jt on the link dual to t using 
the formula 

V(2it + l)Tr \d^^{U^w)- ■ = 5(C7^w ■ ■ -t/^w, 1). 

' L 1 "t J 1 "t 

jt 

If we regard {Ut}t -.tetrahedra as giving a set of transition functions of a principal G-bundle 
over the combinatorial manifold, the condition U„(t) ■ ■ ■ U„(t) = 1 derived here implies that 
the holonomy around the triangle t is trivial. 

Now we can examine the behavior of Zc under the five basic moves of the simplicial 
complex. In the first move, the subdivision of a 4-simplex into five 4-simplexes creates ten 
triangles and ten tetrahedra. Correspondingly we introduce ten spin variables ji, jio on 
the triangles and ten group elements Ui, ...,t/io on the tetrahedra, which give transition 
function of a G-bundle in the interior of the original 4-simplex. By summing over jVs, 
we obtain ten 5-functions, and they restrict the G-bundle in the interior to be fiat. Then 
the ten group elements can be set to 1 by redefining the group elements on the original 
4-simplex. The conditions imposed by these ten 5-functions are not independent and four 
of them are redundant, giving rise to a factor 5(1, 1)^. Thus Zq transforms into 5(1, V)'^Zc 
under the first move. 

If G is a finite group, 5(1, 1) is simply a number of elements in G. But it is divergent 
for a continuous group like SU{2). To regularize this divergence, one may cut off the 
summation as ji < L for some large L as in |0] . The divergent factor is removed from the 



partition function by multiplying 5(1, 1) for each 0-simplex and dividing by 5(1, 1) for each 
1-simplex (note that the numbers of the 0- and 1-simplexes increases by one and five in 
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the first move), and we send L — > oo after tliis. In tliis regularization, tlie invariance under 
tlie first move does not fiold at a finite value of L and it is recovered only after taking 
the limit L —>■ oo. Alternatively one may use the g-deformed version of the model defined 
below. When g is a root of unity, q = e27ri/(fc+2)^ ^^le sum over j's is naturally restricted 
as J < k/2 and the g-analogue of 5(1, 1) is given by 

k + 2 



5(1,1) 



In the so-called BF model ||17|| ||18|| which we will identify as a continuum limit of our 
lattice model, this divergent factor would correspond to a part of its gauge volume. 

In the second move, four triangles are added in the interior of two 4-simplexes and 
summing over spins on them generates four 5-functions, three of which are independent. 
Thus Zc transforms into 5(1, l)Zc under the second move. In this move, the number of the 
0-simplexes does not change while the 1-simplexes increase by one. Thus the divergence is 
removed in the same way as in the first move. On the other hand, no redundant 5-function 
appears in the third move, and Zc remains invariant under this move. The fourth and 
the fifth moves are the inverses of the first and the second moves, and Zc transforms into 
5(1, l)~^Zc and 5(1, l)~^Zc respectively. Combining these result, the dependence of Zc 
on C is expressed as 

Zc = 5(1,1)-^o(^)+^^(^)xg(C), (12) 

where No{C) and Ni{C) are numbers of 0- and 1-simplexes in C, and xciC) depends only 
on the combinatorial class of C. 

It is straightforward to define the g-deformed version of the model, by replacing the 
15j-symbols and the 6j-symbols by their g- analogues. Since the invariance under the 
Alexander moves is essentially due to the Biedenharn-EUiot identity and the orthonormality 
of the 6j-symbols, which are also satisfied by their g-analogues j^], it is likely that the q- 
analogue model is topological. The detailed analysis will be reported elsewhere. 

The expression ( pT|) suggests that, for a compact group G, Zc is related to the partition 
function of the BF model which is defined for an oriented manifold M as 



^BF 



[dB,dA]exp{i [ {B , dA + [A, A])) , 



where A is a connection one-form for a G-bundle on the manifold M, S is a ^-valued two- 
form (^ is a Lie algebra of G) and (, ) is the invariant bilinear form on ^. If we integrate 
over B first, we obtain 

Zbf = j[dA]S{dA+[A,A]), (13) 



and Zbf is given by an integral over flat connections. As we have seen in the above, the 
set of group elements {UT}T:tetrahedra in (01) deflne a principal G-bundle on the simplicial 
manifold, and the sum over j's restricts the G-bundle to be flat. Due to the topological 
nature of the lattice model, we can take the triangulation to be inflnitely flne. When C 
is a smooth triangulation of M, the integral over Ut^s in (|Tl|) is then identifled with the 
integral over A in (|T3| ) in this limit in the sense of the lattice gauge theory. One may also 
regard the spin-j associated to each 2-simplex as a discretized version of the two-form B. 
The action for the BF model has the following gauge symmetries, 

6A^DX , SB = [B, A] + Duo, (14) 

where A and lj are ^-valued 0- and 1-forms, and D is the covariant derivative with respect 
to the connection A. The Wilson operator associated to a closed loop 7 is invariant under 
these gauge transformations and thus gives a physical observable in the theory. One may 
also try an operator deflned by an integral over a closed surface E in M as exp(i /^-.(a, B)) 
for some a in Q. Insertion of this operator generates a holonomy e*" around the surface 
E. This is because the fleld strength F becomes non- vanishing on E upon integration over 
B. In the standard Yang-Mills theory, such an operator has been considered by 't Hooft 
PQI , so we call it the 't Hooft operator here. If we take e*" to be in the center of G, the 
operator becomes invariant under the gauge transformations (0). 



In order to maintain the gauge symmetries without restricting e*" to be in the center, 
we consider the following modiflcation to the 't Hooft operator B. Let us choose a normal 
vector fleld on E and displace E slightly in the direction of the vector fleld to obtain another 
surface E'. The trajectory swept out by E then deflnes a three-dimensional submanifold 
V£ in M bounded by E and E'. This procedure is similar to the framing of the Wilson-line 
in the three-dimensional Chern-Simons theory [^. The integration of the B fleld in Vj: 



restricts A to be flat there, and this induces a flat connection on E as the boundary value 
of A. Especially when E is homeomorphic to S"^, the flat gauge fleld is trivial and the 
parallel transport 0(x, y) = Pexp(z A) is deflned globally on E. Thus, in this case, we 
can deflne a gauge-invariant operator by 

T«(E) = jdUexp{i j (fx{B{x),VL{x,y)UaU-^VL{y,x))). (15) 



^ An operator associated to a surface has also been considered in [24|, but it differs from the 
one constructed here. 
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Because of the [/-integral, this definition of Tq,(E) is independent of the choice of a point 
y on E. It is easy to show that this is invariant under both of the gauge transformations 
in the above. In the BF model, the diffeomorphism invariance is the consequence of these 
gauge symmetries. Thus the expectation value of Tq,(E) is invariant under homotopy moves 
of E. If E is homeomorphic to a torus or a higher-genus surface, the fiat gauge field is not 
necessarily trivial and this construction would need to be modified. In the following, we 
restrict ourselves to the case when E is homeomorphic to S'^. 

Now let us discuss how to realize Tq,(E) in the lattice model. In the expression (|ll|), 
each group element Ut is associated to a tetrahedron T in the simplicial complex C, and 
each 2-simplex carries a spin jf. For each 2-simplex t, the group elements in the tetrahedra 
sharing t are multiplied together in the spin-jt representation D^^ and the trace of the 
product is taken. Therefore the holonomy e*" around a given 2-simplex can be imposed 
by inserting e*" in the trace. In the combinatorial manifold, the surface E is represented 
by a collection of 2-simplexes E. Thus we would like to insert e'" in the traces around all 
the 2-simplexes on E in such a way that the topological invariance is maintained. This can 
be done, by discretizing the expression (p!5[), as follows. 

Let us choose a subcomplex homeomorphic to 5"^ x [D"^: two-dimensional disc) 
in C such that the surface E ~ S"^ is in the image of its boundary (~ S'^ x S^). We first 
sum over the spin-j's on the 2-simplexes in the interior of Vs. The G-bundle in the interior 
becomes flat, and the group elements in the tetrahedra in the interior can be set to 1 by 
redefining the group elements on the boundary of V^. We then insert e*^" on the links dual 
to the 2-simplexes on E. Each 2-simplexes on E is shared by a pair of tetrahedra in dVj:, 
and the insertion of e*" is made on the links connecting these pairs. 

After inserting e*" on E, we integrate over the group elements in the tetrahedra in C. 
In the exterior of Vz, the integral gives the product of the 15j-symbols and the 6j-symbols 
as in ([To|) . Therefore the product becomes invariant under the Alexander moves in the 
exterior of Vs after summing over spin-j's. On the other hand, the integral over the group 
elements on the boundary dVj: generates a product of the 3j-symbols corresponding to the 
(p^-graph on dV^. The product is summed over the magnetic angular momenta on the links 
with insertions of e*" on the links dual to the 2-simplexes on E. This defines a function 
of spin-j's on the ^'^-graph on dVj: ~ S"^ x S"^. By construction, the function depends on 
the triangulation on S'^ x S^, but its dependence is cancelled by the contribution from the 
exterior of Vj: after summing over the spin-j's on dV^. The topological invariance is then 
maintained. 
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In [|T0[ [|T3[, the notion of physical states has been introduced in the topological lattice 
model in three dimensions. In the three-dimensional model, states on a two-dimensional 
closed surface are functions of spin-j's on a (^'^-graph dual to a triangulation of the surface, 



and they are called physical when they transform in a definite way specified in [|T3[ under 
the Alexander moves of the triangulation. It is possible to introduce a similar notion in 
the four-dimensional model here. States in the four-dimensional model are functions of 
spin-j's on a (p^-graph on a three-dimensional closed manifold, and they are called physical 
when they transform in such a way that the insertion of the state on the three-manifold 
embedded in four dimensions does not spoil the topological invariance. We may regard 
the 't Hooft operator as creating a physical state on S"^ x S^. 

For the 't Hooft operator, S"^ x is the boundary of Vy, S"^ x D"^ , and the S'^-part 
is contractable in M. One may also consider a complementary situation when the iS^-part 
is contractable while the S'^-part may link with other operators. This gives the Wilson 
operator for a closed loop 7 which is the image of . To be more explicit, we consider a 
subcomplex homeomorphic to x [D^: three-dimensional ball) in C such that the 
loop 7 travels through tetrahedra Ti,...,T^ in the interior of W-y. We then make a choice of 
a spin j and insert the discretized version of the Wilson operator tr[D^ {Uti) ■ ■ ■ D^{Ut„)] 
in the integrand of (pT|). By integrating out the group elements in the tetrahedra in Wj 
and by summing over the spins on the links in the interior of W^, we obtain a function of 
the spins on the boundary of W^. It is easy to check that the function defines a physical 
state on 5"^ X 5"^. 

We have constructed the 't Hooft operator associated to a surface E and the Wilson 
operator associated to a loop 7. They both create physical states on x S^, so their 
presence does not spoil the topological invariance of the lattice model. Especially their 
correlation functions are invariant under ambient isotopy moves of E and 7, and thus give 
invariants of knots and links. It should also be possible to define the 't Hooft operator 
and the Wilson operator in the g-deformed version of the lattice model. Detailed study of 
these invariants is left for future investigations. 
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